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As is well known, the operator product expansion (OPE) and QCD sum rule method 
can be used in QCD for problems characterized by short distances much smaller than the 
QCD scale, / ^ A ~ 0.5 fm. Under no circumstances can these methods be exploited to 
treat long distance problems, / > 1 fm. As was shown in our letter, meson (and baryon) 
mass shifts Am in nuclear matter are determined by long distances. This follows from the 
basic relation of our letter 

Am{E) = -27r£Re/(^) (1) 

in which the mass shift Am is expressed in terms of the real part of the meson-nucleon 
forward scattering amplitude. Re/, and where p is the nuclear density. Wc showed that 
for p- mesons |Re/| > 1 fm. The same is true for 7r-mesons at £J > 170 MeV. If it were 
possible to determine Am in the short distance expansion, then Re/ > 1 fm could also 
be obtained in this way, which is obviously not true. It must be stressed, that Eq. (1) is 
valid starting from very low energies. The only restriction comes from the requirement that 
the particle wavelength be less than the intcrnuclcon distance d. For p-mcsons and normal 
nuclear density this restriction gives E^^"' > 10 MeV. In our consideration of the p-meson 
mass shift we confined ourselves to high energies, ^1.2 GeV, only because we could not 
find Re/pjv at lower energies using the vector dominance model. (This restriction is absent 
in the pion case where we use experimental data on Re/^Ar). Since the region E^"' < 1 GeV 
is the resonance region in pA^-scattering, we do not believe that, in this region and near the 
threshold, Re/pjv is much less than 1 fm. However, should it somehow happen to be the case 
(with short distance methods becoming apphcable), then we would get |Amp| < 10 MeV, 
which still would not agree with the sum rule results of Hatsuda and Lee. 

Let us remark that the two versions of OPE mentioned by Hatsuda and Lee, the light- 
cone (LC) expansion and the short distance (SD) expansion, are not in fact two different 
options. The second is a limiting case of the first. In the space-time picture both correspond 
to expansion near the light cone, the second being done in the vicinity of its tip. Which 
version to choose is not in our hands - it is dictated by the physics of the problem under 
consideration. In the case of a two-body problem (and the meson mass shift in nuclear 
matter is a two- or even many-body problem) the longitudinal distances along the light cone 
are of order I ~ 1/mx, where m is the nucleon mass and x — Q^/2i/ is the Bjorken variable. 
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In case of a p-meson at rest ~ m^, v ~ mnip, and I is rather large, I ~ 2/mp ~ 0.5 fm. 
Therefore, the SD expansion is not apphcable here. It must also be mentioned, that in using 
fixed-q dispersion relations, like the one in Eq. (3) of the Comment, a serious problem arises 
in the phenomenological part of the sum rule. One must separate the contributions of states 
in the meson channel which are of interest from the contributions of baryonic resonances in 
s- and, especialy, li-channels. The results of such separation are strongly model-dependent. 

The statement made in the Comment that we overlooked in our letter the mean-field 
(MF) description of low-q mesons is due to a misunderstanding which perhaps stems from 
the fact that we did not use the term "mean field theory" explicitly. (However, this term 
was used in Refs. [16,17] we quoted in our derivation). In fact, our basic formula, Eq. 
(1), is a direct consequence of the MF theory. The derivation of this formula in our letter 
used MF concepts formulated in coordinate space and not in, the more commonly used, 
momentum space. It can also be derived from Eq. (1) of the Comment by considering the 
pole contribution to the l.h.s. of that equation. In our case, when the meson momentum 
is much larger than the nucleon momenta in matter, the r.h.s. of Eq. (1) of the Comment 
factorizes and reduces to the nucleon density appearing in our Eq. (1). Of course, at a low 
meson momentum one should integrate over the nucleon momenta taking into account that 
R,ef{E) is a function of the total energy in the cm. frame. Positions and properties of 
resonances in meson-nucleon scattering are very important in such a calculation, although 
in other aspects it is trivial. At very low meson momenta, |q|~^ ^ d, the averaging over the 
Fermi sea does not make sense anymore, since interactions with more than one nucleon and 
screening effects become important, and the MF approximation in its traditional form (linear 
in the Fermi distribution) breaks down. These effects are not accounted for in calculations 
based on the SD expansion. 
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